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ABSTRACT 

As is well known, symmetry plays an important role in the theoretical physics. In particular, the 
well-known Noether symmetry is an useful tool to select models motivated at a fundamental level, 
and find the exact solution to the given Lagrangian. In the present work, we try to consider Noether 
symmetry in f(T) theory. At first, we briefly discuss the Lagrangian formalism of f(T) theory. In 
particular, the point-like Lagrangian is explicitly constructed. Based on this Lagrangian, the explicit 
form of f(T) theory and the corresponding exact solution are found by requiring Noether symmetry. 
In the resulting f(T) = fj,T n theory, the universe experiences a power-law expansion a(t) ~ t 2 ™/ 3 . 
Furthermore, we consider the physical quantities corresponding to the exact solution, and find that 
qq ■ if n > 3/2 the expansion of our universe can be accelerated without invoking dark energy. Also, 

we test the exact solution of this f(T) theory with the latest Union2 Type la Supernovae (SNIa) 
dataset which consists of 557 SNIa, and find that it can be well consistent with the observational 
data in fact. 
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I. INTRODUCTION 

The current accelerated expansion of our universe [l[ has been one of the most active fields in modern 
cosmology. As is well known, it could be due to an unknown energy component (dark energy) or a 
modification to general relativity (modified gravity) 0, Q • Today, modified gravity theory has been a 
competitive alternative to the familiar dark energy scenario. 

In analogy to the well-known f(R) theory, recently a new modified gravity theory, namely the so-called 
/ (T) theory, has been proposed to drive the current accelerated expansion without invoking dark energy. 
It is a generalized version of the so-called teleparallel gravity originally proposed by Einstein 0, Q ■ In 
teleparallel gravity, the Weitzenbock connection is used, rather than the Levi-Civita connection which is 
used in general relativity. Following || , here we briefly review the key points of teleparallel gravity 
and f(T) theory. In this work, we consider a spatially flat Friedmann- Robertson- Walker (FRW) universe 
whose spacetime is described by 

ds 2 = -dt 2 + a 2 (t)dx 2 , (I) 

where a is the scale factor. The orthonormal tetrad components e^x^) relate to the metric through 

g^u = m^el , (2) 

where Latin i, j are indices running over 0, I, 2, 3 for the tangent space of the manifold, and Greek 
p, v are the coordinate indices on the manifold, also running over 0, I, 2, 3. In teleparallel gravity, the 
gravitational action is 

S T = J d i x\e\T, (3) 

where |e| = det (e^) = \J—g , and for convenience we use the units 16irG = fo = c = 1 throughout. The 
torsion scalar T is defined by 



T = S i»> T p ( 4 ) 



where 



TV = ~< (^4 - d v ei) , (5) 
K»\ = ~ (T"% - T v % - T/») , (6) 

Sr = \ {K^ p + S;T ev e - S;T%) . (7) 

For a spatially flat FRW universe, from Eqs. (HJ) and ([!]), one has 

T = -QH 2 , (8) 

where H = a/ a is the Hubble parameter, and a dot denotes a derivative with respect to cosmic time t. In 
analogy to the well-known f(R) theory, one can replace T in the gravitational action ([3]) by any function 
/(T), and then obtain the so-called f(T) theory. In f(T) theory, the modified Friedmann equation and 
Raychaudhuri equation are given by [5H8J 

l2H 2 f T + f = p, (9) 
48H 2 f TT H - h (l2# 2 + 4if) -f=p, (10) 

where a subscript T denotes a derivative with respect to T, and p, p are the total energy density and 
pressure, respectively. In an universe with only pressureless matter, obviously we have p = p rn = and 
P = Pm = PmO a ~ 3 , where the subscript "0" indicates the present value of the corresponding quantity, 
and we have set ao = I. It is well known that when f(T) = T the familiar general relativity can be 
completely recovered. 



3 



In fact, f(T) theory was firstly used to drive inflation by Ferraro and Fiorini Q. Later, Bengochea and 
Ferraro [5[, as well as Linder [6j, proposed to use f(T) theory to drive the current accelerated expansion 
of our universe without invoking dark e nerg y. Very soon, f(T) theory attracted much attention in the 
community We refer to e.g. @, 13, E3-EH E5 for relevant works. 

So far, the specified forms of function /(T) in the literature are written by hand. There is no natural 
guidance from fundamental physics on the form of /(T) in fact. In the present work, we try to address this 
issue. As is well known, symmetry plays an important role in the theoretical physics. In particular, the 
well-known Noether symmetry is an useful tool to select models motivated at a fundamental level. In the 
literature, Noether symmetry has been extensively used in scalar field cosmology [lfl non-mini mally 
coupled cosmology fl3. f(R) theory fl^-pll . scalar-tensor theory (22|, higher order gravity theory 23 1, 
multiple scalar fields 1241 . vector field [25|], fermion field [26|], tachyon field (27|, non-flat cosmology |28l |. 



quantum cosmology |29|, Bianchi universe l30|, Brans-Dicke theory [31|, dilaton [32], induced gravity 
theory (33[, gravity with variable G and A |34{ . Gauss-Bonnet gravitvj35| . and so on. Now, we try to 
consider Noether symmetry in f(T) theory in the present work. 

This paper is organized as follows. In Sec. HH we briefly discuss the Lagrangian formalism of f(T) 
theory. In particular, the point-like Lagrangian is explicitly constructed. In Sec. lIIIl we consider Noether 
symmetry in f(T) theory. The explicit form of f(T) theory and the corresponding exact solution are 
found by requiring Noether symmetry. In Sec. lIVl we further discuss the physical quantities corresponding 
to the exact solution. In Sec.|Vl we test the exact solution of f(T) theory found in Sec. lIIIl with the latest 
Union2 Type la Supernovae (SNIa) dataset which consists of 557 SNIa. Finally, some brief concluding 
remarks are given in Sec. I VII 



II. LAGRANGIAN FORMALISM OF f(T) THEORY 

In the study of Noether symmetry, the point-like Lagrangian plays an important role. In this section, 
we discuss the Lagrangian formalism of f(T) theory. As mentioned above, the relevant action of f(T) 
theory is given by 



S = 



J d A x\e\f{T)+S„ 



(11) 



where S m is the action of pressureless matter minimally coupled with gravity, and we assume that the 
radiation can be ignored. Following e.g. [2(| [2l|, to derive the cosmological equations in the FRW metric, 
one can define a canonical Lagrangian C = £(a,d,T, T), whereas Q = {a,T} is the configuration space, 
and TQ = {a, a, T, T} is the related tangent bundle on which £ is defined. The scale factor a(t) and the 
torsion scalar T(t) are taken as independent dynamical variables. One can use the method of Lagrange 
mutipliers to set T as a constraint of the dynamics (nb. Eq. (|5|)). Selecting the suitable Lagrange mutiplier 
and integrating by parts, the Lagrangian C becomes canonical [20l |21|. In our case, we have 



5 = 2^ 



dta 6 



f(T)-X(T + 6^)-^f 



where A is a Lagrange mutiplier. The variation with respect to T of this action gives 

A = h ■ 

Therefore, the action (fT2"|) can be rewritten as 



S = 2it z / dt a 



f{T)-h [T + 6 



PmO 
n3 



(12) 



(13) 



(14) 



and then the point-like Lagrangian reads (up to a constant factor 2ir 2 



C(a, a, T, t) = a 3 (/ - f T T) - Qf T aa 2 - p m0 . 



(15) 
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As is well known, for a dynamical system, the Euler-Lagrange equation is 

m- a i=°> 

where qi are the generalized coordinates of the configuration space Q, and in our case qi = a and T. 
Substituting Eq. (TT5|) into the Euler-Lagrange equation (|16p . we obtain 

a 3 f TT ( T + 6 J) =0, (17) 



/ - f T T + 2f T H z + 4 I / T - + ff/rrT 1 = . (18) 
If /tt 7^ 0, from Eq. (TT71) . it is easy to find that 

. 2 

T=-6% = -6H 2 , (19) 
or 

i.e., the relation ([8]) is recovered. Generally, this is the Euler constraint of the dynamics. Substituting 
Eq. (HU) into Eq. (JTHJ) and using d/a = H 2 + H, we get 

A8H 2 f TT H - 4/ T (SH 2 + H) - f = , (20) 

i.e., the modified Raychaudhuri equation (|10[) is recovered (note that p = p rn = 0). On the other hand, 
it is also well known that the total energy (Hamiltonian) corresponding to Lagrangian C is given by 

n = T,7T^- C - ( 21 ) 
Substituting Eq. $TQ) into Eq. (HTJ), we have 

H(a, a, T, t) = a 3 ( -6f T ^ - f + f T T + ^) . (22) 



Considering the total energy H = (Hamiltonian constraint) |16l. |20L \2l\ and using Eq. (flU)) . we get 

12H 2 f T + f = ^ , (23) 

i.e., the modified Friedmann equation Q is also recovered (note that p = p rn = p m o/a 3 ). So far, we have 
shown that the point-like Lagrangian given in Eq. (|15p can yield all the correct equations of motion, and 
hence it is the desired one. 



III. NOETHER SYMMETRY IN f(T) THEORY 



As is well known (see e.g. [1 61 - 4351 ] ^ . Noether symmetry is an useful tool to select models motivated at a 
fundamental level, and find the exact solution to the given Lagrangian. In this section, we try to consider 
Noether symmetry in f (T) theory. 



Following e.g. [20L l21j . the generator of Noether symmetry is a vector 



d d d ■ d 
X = ^ + ^T + ^ + ^?< (24) 



where a = a(a,T) and f3 = /3(a,T) are both functions of the generalized coordinates a and T. Noether 
symmetry exists if the equation 



(25) 
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has solution, where L^C is the Lie derivative of the Lagrangian C with respect to the vector X. Of course, 
according to the well-known Noether theorem, there will be a constant of motion (Noether charge) (20ll2lj. 
namely 

_ x - dC dC n dC . . 

Qa = } ^i— = a — + fj—r- = const. (26) 
^ dqi da dT 

Note that the meaning of Ly^L = is that C is constant along the flow (possibly a local flow) generated 
by X, namely Eq. ((25)) is identically verified all over TQ [l6[. Its explicit evaluation gives an expression 
of second degree in a and T, whose coefficients are functions of a and T only. Therefore, they should be 
zero separately [lg, |2(J, |2l| . 

In our case, substituting Eq. (TTSjl into Eq. (|2"5"]) and using the relations a = (da /da) a + (da/dT)T, 
P = (df3 /da) a + (dp/dT) f, we obtain 

3aa 2 (f - f T T) - f3a 3 f TT T - 6a 2 (af T + Mtt + 2a fr^j ~ = . (27) 

As mentioned above, requiring the coefficients of a 2 , T 2 and dT in Eq. (|2"T)) to be zero, we find that 

4f =o - (28) 

da 

af T + f3af T T + 2af T -^=0, (29) 

3aa 2 (f - f T T) - (3a 3 f TT T = . (30) 

In particular, the constraint (|3"0)) is sometimes called Noether condition [20(. The corresponding constant 
of motion (Noether charge) given in Eq. (|2l)]) reads 

Qo = —\2afTaa — const. (31) 

A solution of Eqs. (|28p. ([2^)1 and (T5TTT) exists if explicit forms of a and /3 are found; and if at least one 
of them is different from zero, Noether symmetry exists [20(. Obviously, from Eq. (|28|). it is easy to see 
that a is independent of T, and hence it is a function of a only, i.e., a = a(a). On the other hand, from 
Eq. ([3"0"]l. we have 

Paf TT T = 3a (/ - f T T) . (32) 
Multiplying T in both sides of Eq. (|2"9")l , and then substituting Eq. (|3"2")l into it, we obtain 

f T T (%a^ - 2a^j +3af = 0. (33) 

Fortunately, one can perform a separation of variables, and recast Eq. (|33[) as 

l _ada = Jf 
a da 2jtT 

Since its left-hand side is a function of a only and its right-hand side is a function of T only, they must 
be equal to a same constant in order to ensure that Eq. always holds. For convenience, we let this 
constant be 3/(2n), and then Eq. (134[) can be separated into two ordinary differential equations, i.e., 

nf = f T T, (35) 

(36) 

a da 2n 

It is easy to find the solutions of these two ordinary differential equations, namely 

f(T) = fiT n , (37) 
a(a) = a a 1 - 3 '^ , (38) 
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where p, and «o are integral constants. Obviously, f(T) and a{a) are both power-law forms. Substituting 
Eqs. dSTJ) and {33) into Eq. K32J), we find that 

^ (aj r) = -^ a - 3 /( 2n )T. (39) 

n 

So far, we have found the explicit non-zero solutions of f{T), a and /?. Therefore, Noether symmetry 
exists. Finally, we try to find out the exact solution of a(t) for this type of f(T). Substituting Eqs. (|3"T)) . 
and (1191) into Eq. (|31[) . we can obtain an ordinary differential equation of a(t), namely 



a Cl a = c 2 , (40) 

where 



3 1 



Qo 



-i l/(2n-l) 

-12a Ain(-6)"- 1 _ 
It is easy to find that the solution of Eq. P0"j) is given by 



(41) 



0(t) = -(1 + CX)(C 3 - C2 t) 1 /d+«=i) = (-1)1+W3 . 3 {c2t _ C3 )2n/3 (42) 

where C3 is an integral constant. Obviously, in the late time \c2t\ ^> | C3 1 and the universe experiences a 
power-law expansion. In fact, we can make it clearer. Requiring a(t = 0) = 0, it is easy to see that the 
integral constant C3 is zero in fact. So, we have 



a(t) ~ i 2n/3 , (43) 
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where its prefactor (— l) 1+2n / 3 . __ c% is not important. Note that n > is required to ensure that the 
universe is expanding. 



IV. PHYSICAL QUANTITIES CORRESPONDING TO THE EXACT SOLUTION 

Here, we further consider the physical quantities corresponding to the exact solution found in the 
previous section. Firstly, from Eq. (1431) . we find that the Hubble parameter is 

H=- = —t-\ (44) 
a 6 



and the deceleration parameter is given by 



aH 2 2n 



£-1- ( 45 ) 



When n > 3/2, the expansion of our universe can be accelerated (note that n < is not acceptable since 
the universe contracts in this case). Secondly, from Eqs. ([9]) and (fT0|) . one can find that the effective dark 
energy density and pressure from torsion are given by [5|-|8( 

Pde = 6H 2 -f- 12H 2 f T , (46) 
Pde = -Pde - 4 (l2H 2 f TT -f T + l)H. (47) 

Substituting Eqs. (|gi)). (TH?|) and (J37J) into Eqs. (gSJl and P7|). we can find that the equation-of-state 
parameter (EoS) of the effective dark energy from torsion is given by 

= Pde_ ^ &n(n - 1) 3" 

Wde ~ p de 8n 2 • 3" - 3/i(8" - n ■ 2 1 + 3n )(-n 2 ) n t 2 ^) ' 1 ' 
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Using the well-known relation between the total EoS Wtot and the deceleration parameter q (see e.g. [36j), 
it is easy to see that 

Wtot =P^ = l(2q-l) = --l. (49) 
Ptot o n 

Again, when n > 3/2 we have Wtot < — 1/3 and then the expansion of our universe can be accelerated 
(note that n < is not acceptable since the universe contracts in this case). It is worth noting that if 
n > 1, from Eq. (|48[) we have Wde 1M — 1 = Wtot hi the late time t — > oo. This is not surprising. Note 
that the well-known relation (see e.g. J36|) 

Wtot = MdeWde + SlmW m , (50) 

where the EoS of pressureless matter w m — 0, and fii are the fractional energy densities of the effective 
dark energy and pressureless matter. Although u> tot is constant, Wde and £ld e are dependent on time t. 
If Wde < w m = 0, it is inevitable that Qd e — > 1 and then Wde w to t in the late time t —¥ oo. Finally, 
we turn to the fractional energy densities. Using Eqs. (|4"8"| . (|4"9"|) . ([ST))) and w m = 0, we find that the 
fractional density of the effective dark energy from torsion is given by 

w tot _ 3^(8" -n-2 1 + 3 ")(-n 2 )"t 2 ( 1 -") 

S2de -^7 _1 8^^ ' ^ 

and then the fractional density of pressureless matter J7 m = 1 — Qd e is ready, namely 

3/^(8" -n-2 1+3Tl )(-n 2 ) n i 2 ( 1 -™' , , 

^ = 8^ ■ (52) 

It is worth noting that f2 m coming from Eq. (|52[) contains only two parameters n and /Lt. On the other 
hand, by definition Q m oc p m /H 2 and /5 m = p m o/a 3 , while the prefactor in a(t) contains C2 which is given 
in Eq. (|4ip . we see that VL m contains five parameters p m o, n, fi, ao and Qo in this way. Requiring the 
equality between these two fi m , it is easy to see that at least one of these five parameters p m o, n, (i, ao 
and Qo is not independent. 



V. COSMOLOGICAL TEST 



As is well known, the current accelerated expansion of our universe was firstly found from the obser- 
vation of distant Type la Supernovae (SNIa) [l[. Here, we would like to test the exact solution of f(T) 
theory found in Sec. Mil with the latest Union2 SNIa dataset [37j which consists of 557 SNIa. 

The data points of the 557 Union2 SNIa compiled in [32| are given in terms of the distance modulus 
p, bs{zi). On the other hand, the theoretical distance modulus is defined as 

Pth(zi) = 51og 10 L» L (z. i ) + , (53) 

where z = 1/a — 1 is redshift, po = 42.38 — 51og 10 /i and h is the Hubble constant Ho in units of 
100 km/s/Mpc, whereas 

D ^»- (1+ "jCifey (54) 

in which E = H/Hq, and p denotes the model parameters. Correspondingly, the \ 2 from the 557 Union2 
SNIa is given by 



where a is the corresponding \a error. The parameter \io is a nuisance parameter but it is independent of 
the data points. One can perform an uniform marginalization over (1q. However, there is an alternative 
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way. Following |38l 1391 ] , the minimization with respect to can be made by expanding the xfi of Eq. (|55[) 
with respect to /^o as 



where 



X 2 Jp)=A-2 l x B + f 4C : 



Mp) = E 



(56) 



*(p) = E 



fJ-obsjzj) - Hth{zi\ Mo = 0,p) 



C = 



,(*) 



Eq. (j56|) has a minimum for /ip = B/C at 



^(P) = i(p)-^ 



(57) 



Since xfi min = x 2 m in obviously (up to a constant), we can instead minimize xft which is independent 
of hq. As is well known, the best-fit model parameters arc determined by minimizing \ 2 — xt- As 



Xmin — 1-0) 2-3 and 3.53 for n p 



1. 



in [3j||40j], the 68.3% confidence level is determined by Ax" = x~ 
2 and 3, respectively, where n p is the number of free model parameters. Similarly, the 95.4% confidence 
level is determined by Ax 2 = X 2 ~ Xmin — 4-0, 6.17 and 8.02 for n p = 1, 2 and 3, respectively. Note that 
the corresponding h can be determined by (1q = B/C for the best-fit parameters. 





FIG. 1: The X 2 and likelihood £_, 2 oc e x 7/2 as functions of parameter 



In our case, from Eqs. (|4"3"]) and (|44D . we have H oc a 3 /( 2 "), and then 

E= — = a - 3 /( 2 ») = (1 + zf/^ . 
Ho 



(58) 



We present the corresponding x 2 and likelihood C x 2 oc e x I" 1 as functions of parameter n in Fig. [T] The 
best fit has Xmin = 549.222, and the best-fit parameter is 

n = 2.417±ai3i (with la uncertainty) IqJso (with 2a uncertainty). (59) 
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The corresponding h — 0.692 for the best fit. Obviously, n > 3/2 as required by the accelerated expansion 
of our universe. In Fig. [21 we show the Hubble diagram for the best fit, comparing with the 557 Union2 
SNIa data points. Obviously, our f(T) theory with solution (PfBI can be well consistent with the latest 
557 Union2 SNIa dataset. 

46 
44 
42 
40 

a. 

38 
36 
34 



0.2 0.4 0.6 0.8 1 1.2 1.4 

z 

FIG. 2: The Hubble diagram for the best fit (red solid line), comparing with the 557 Union2 SNIa data points 
(black diamonds). See the text for details. 

VI. CONCLUDING REMARKS 

As is well known, symmetry plays an important role in the theoretical physics. In particular, the well- 
known Noether symmetry is an useful tool to select models motivated at a fundamental level, and find the 
exact solution to the given Lagrangian. In the present work, we try to consider Noether symmetry in f(T) 
theory. At first, we briefly discuss the Lagrangian formalism of /(T) theory. In particular, the point-like 
Lagrangian is explicitly constructed. Based on this Lagrangian, the explicit form of f(T) theory and the 
corresponding exact solution are found by requiring Noether symmetry. In the resulting f(T) = /iT" 
theory, the universe experiences a powerdaw expansion a(t) ~ t 2 ™/ 3 . Furthermore, we consider the 
physical quantities corresponding to the exact solution, and find that if n > 3/2 the expansion of our 
universe can be accelerated without invoking dark energy. Also, we test the exact solution of this f{T) 
theory with the latest Union2 SNIa dataset which consists of 557 SNIa, and find that it can be well 
consistent with the observational data in fact. 

Some remarks are in order. Firstly, although there are many works concerning Noether symmetry in 
the literature, we still find something new by considering Noether symmetry in f(T) theory. For instance, 
we take f(R) theory for comparison. It is well known that the equations of motion in f(R) theory are 
4th order while they are 2nd order in /(T) theory; R contains a and a, but T contains only a. Therefore, 
the corresponding point-like Lagrangian of f(R) theory contains R, but the one of f{T) theory does 
not contain T. Note that in the Ly^L = equation, the Euler-Lagrange equation and the total energy 
(Hamiltonian) corresponding to Lagrangian, dC/dqi plays an important role. The absence of T in the 
point-like Lagrangian of f(T) theory makes difference and brings something new. Secondly, unlike other 
theories (e.g., f(R) theory, which has been studied for more than ten years), f(T) theory became an active 
field just from 2010. Due to its relatively short history, many aspects of f(T) theory are unexplored in 
fact. For instance, to our knowledge, the reconstruction of f(T) from cosmological observations is absent 
in the literature so far, unlike f(R) theory. We do not know what forms of f(T) are the suitable ones, or 
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what forms of f(T) have a good motivation from fundamental theories and principles. So far, all f(T) 
forms considered in the literature are written by hand. In such a situation, it is important to find a f(T) 
form from some principles. Of course, one might argue that Noether symmetry is not really a physical 
guiding principle. But, why we deny to use Noether symmetry to shed some light on the unknown side 
of f(T) theory? So, we consider that it is better to keep an open mind. Finally, we have shown in Sec. fVl 
that the resulting f(T) theory from Noether symmetry can be well consistent with the observational 
data. We are not playing mathematical games with Noether symmetry. Instead, we have shown that the 
f(T) motivated by Noether symmetry can be tested by realistic cosmological data, and then it can be 
considered as a serious theory of modified gravity and an alternative to dark energy. 
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